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First Observation of Accelerator Muon Antineutrinos in MINOS 
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We report the first direct observation of muon antineutrinos in the MINOS Far Detector in the current muon- 
neutrino dominated beam. The magnetic field of the detector is utiUzed to separate muon neutrinos and 
antineutrinos event-by-event by identifying the charge sign of the muon created in charged-current interactions. 
We present preliminary results on the oscillation parameters as well as limit on the fraction of neutrinos that 
disappear and reappear as antineutrinos. We also discuss the prospect of the measurement when the polarity 
of the magnetic focusing horns will be reversed to create a dedicated muon antineutrino beam. 



1. Introduction 

MINOS is a long-baseline neutrino experiment in 
the NuMI (Neutrinos at the Main Injector) beam line 
at Fermilab that has confirmed the disappearance of 
muon neutrinos and measured the atmospheric oscil- 
lation parameters (sin^(26') and |Am^|) with high pre- 
cision In this paper, we describe a new study that 
utilizes the 7% muon antineutrino component of the 
beam to measure their oscillation parameters directly 
and test exotic models such as CPT violation in the 
neutrino sector Q and transition of f ^ to Q . 

In a quasi-two-neutrino mixing framework, oscil- 
lation implies a survival probability of 

P(z>,^P,) = l-sin2(20)sin2(^l:^I|^), (1) 

where E is the neutrino energy in GeV, L is the dis- 
tance traveled by the neutrino in km, and the antineu- 
trino mixing angle (9) and mass separation (Am^) are 
assumed to be independent of the corresponding neu- 
trino parameters^ (0 and Arn^). In addition to the 
oscillation scenario, we also consider the possibility 
that a fraction, a, of the muon neutrinos that have 
been observed to disappear along their long flight [ll] 
will reappear as antineutrinos. The probability of this 
transition can be parametrized empirically as 

P(^^. ^ p^) =asin\29)sm' |^ 1.27Am^Z, ^| ^2) 

2. Neutrino Beam and Detectors 

MINOS uses the high intensity neutrino beam [3] 
created by 120 GeV protons from Fermilab's Main 
Injector impinged on a graphite target. Secondary 
particles, mainly tt and K mesons, from the target 



^In this framework, i/^^ oscillation is dominated by the mixing 
between the U2 and 1^3 mass eigenstates, thus Am^ ^™23 ~ 
mg — and 6 fa 623- Similarly for the antineutrinos. 



are sign-selected and focused by two toroidal magnetic 
horns toward a 675 m long evacuated pipe where they 
can decay in flight to produce neutrinos. In the cur- 
rent beam conflguration the horns are focusing (de- 
focusing) secondary particles with positive (negative) 
charge, which enhances the fraction of muon neutri- 
nos at the expense of antineutrinos. Most of the 
arises from tt~ (and K~) parents produced upstream 
in the target that travel down the center of the horns 
where they are not deflected by the magnetic field. 
In addition, a significant portion of the originate 
from parents produced downstream from the target in 
secondary interactions in the decay pipe wall and the 
surrounding material. 

The MINOS experiment [| uses two functionally 
identical scintillator tracking/sampling-calorimeter 
detectors located 1 km (Near Detector) and 735 km 
(Far Detector) from the target. This two-detector 
setup significantly reduces the systematic uncertain- 
ties due to neutrino fiux, cross section, and detection 
efficiency. The detector design is optimized to ob- 
serve charged-current (CC) muon (anti)neutrino inter- 
actions which produce a prominent muon track pen- 
etrating several layers of alternating iron and scintil- 
lator planes. Both detectors have a toroidal magnetic 
field with an average field strength of about 1.3 T 
(near) and 1.4 T (far) in the steel. The momentum of 
the muons is measured from their range in the detec- 
tor or alternatively from their curvature in the mag- 
netic field. In addition, positive (negative) muons are 
focused (defocused) by the field allowing to discrim- 
inate between ^"^-00 and D^-CC interactions event- 
by-event. 

Fig. [T] shows the predicted number of and i?^ 
CC interactions in the MINOS Near Detector and the 
contribution of different beam components to the i?^- 
CC spectrum. Due to the focusing/defocusing effects 
of the horns on the parent mesons, the shape of the 
and spectra are significantly different. The higher 
peak energy of the spectrum around 8 GeV makes 
the measurement more sensitive to higher Am^ 
values than the ly^. 
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Figure 1: Top: the relative size of the S'^-CC energy spec- 
trum (red) compared to the u^-CC spectrum (black) in the 
Near Detector. Bottom: the contribution of the different 
beam components to the Near Detector i^f^-CC spectrum. 



3. Data Analysis 



The data used in this analysis corresponds to 3.2 x 
10^*^ protons on target (PoT) collected in the low- 
energy beam configuration^ between 2005-2007. We 
measure the rate of inclusive i^^i-CC interactions as a 
function of reconstructed antineutrino energy in the 
Near Detector and extrapolate the spectrum to the 
Far Detector, i?^ oscillation would produce an energy- 
dependent deficit while t'^ — > transition would 
cause an excess in the Far Detector compared to the 
expectation. 



^The peak energy of the neutrino beam can be tuned by 
adjusting the relative position of the target and magnetic horns 
as well as by changing the horn currents. 



3.1. Event Selection 

We select events with at least one reconstructed 
track, the longest one identified as the muon candi- 
date. The neutrino interaction point (vertex) is re- 
quired to be inside the fiducial volume in order to sup- 
press cosmic-ray and rock muons originating from out- 
side of the detector as well as to contain the hadronic 
shower in the detector. Cosmic-ray background is 
suppressed further in the far detector by requiring 
the angle of the muon momentum to be within 53° 
(cosi? > 0.6) of the neutrino direction and the event 
time to fall within 14 /is window around the time of 
the beam spill. 

The muon candidate is required to have a positive 
charge determined from the direction of bending in 
the magnetic field. As Fig. [2] demonstrates, this sim- 
ple charge-sign selection produces a sample which is 
highly contaminated. The background is composed 
of both neutral-current (NC) events in which one or 
more tracks are found by the reconstruction code and 
high inelasticity v^-CC events in which a low energy 
/i~ track is obscured by the hadronic shower. In addi- 
tion, as the momentum of the muon increases it bends 
less making the charge determination harder and lead- 
ing to an increasing i/f^ contamination at high energy. 



o 



.9? 2 
o 



LU 



Q 



Far Detector, Low Energy Beam 
MINOS Preliminary 
IVIonte Carlo 




Contamination in Positive Sample 
Contamination after Selection 
Selection Efficiency 



—t—r-U**i \ I irrtTTrTTi 



IE 







5 10 15 20 30 40 50 

v. Energy (GeV) 



Figure 2: The efficiency (solid red) and contamination 
(solid blue) of the final z5,,-event selection according to 
Monte-Carlo simulation of the far detector. The contami- 
nation in the sample with a simple selection on the charge 
of the muon candidate from the track fit is also shown 
(dashed blue). 

Three additional selection variables are used in or- 
der to suppress both the NC and mis-identified Vfj,- 
CC background. The first variable is a likelihood- 
based selector that combines three probability dis- 
tribution functions related to the event topology 
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event length, fraction of total event signal produced 
by the muon candidate, and average signal per plane 
induced by the muon candidate. As Fig. [3] demon- 
strates, this variable is very effective to discriminate 
CC events with a muon track from NC events with 
a diffuse hadronic shower as well as the those j^^-CC 
events that pass the other requirements. The other 
two variables are designed to improve the charge-sign 
determination of the muon candidate. The first vari- 
able is the significance of the measured curvature from 
the fit, and the other is the angle defined in the plain 
transverse to the beam between the direction of the 
last muon hit with respect to the projected hit without 
magnetic field and the line going through the center 
of the magnet coil and the interaction vertex Q ■ 
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The selection was optimized for maximum efficiency 
X purity below f GeV neutrino energy in order to in- 
crease the sensitivity for CPT-conserving oscillation. 
Fig. 121 shows the relative efficiency of the full selection 
and the remaining contamination as a function of 
energy in the Far Detector. According to the MC sim- 
ulation, the overall selection efficiency is 83% with a 
purity of 97% in the Far Detector and f^j-CC events 
are suppressed by a factor of 1.3 x 10"'^. 



3.2. Near-to-Far Extrapolation 

Even in the absence of oscillation or transition, the 
shape of the (anti)neutrino spectra at the Near and 
Far Detector are not identical due the different solid- 
angle coverage of the two detectors and the fact that 
the Near Detector sees an extended source of neu- 
trinos while the Far Detector essentially sees a point 
source. The beam-line geometry and the meson de- 
cay kinematics are encapsulated in a beam-transfer 
matrix that relates the Near Detector energy spec- 
trum to the Far Detector spectrum Using this 
matrix the measured Near Detector spectrum can be 
extrapolated to get a Far Detector prediction, which 
is less sensitive to uncertainties in the neutrino flux 
calculation. Separate transfer matrices are used for 
neutrinos and antineutrinos to predict the signal and 
background simultaneously Q. Simulation of the de- 
tector response is used to correct for energy resolution 
as well as differences in selection efficiency and con- 
tamination between the two detectors. 

In order to improve the extrapolation further, the 
Near Detector data are used to constrain the hadron 
production off the target, which has the largest contri- 
bution to the uncertainty in the neutrino-flux calcula- 
tion. The production of 7r+ and is parametrized 
as a function of their initial transverse (pt) and lon- 
gitudinal (pz) momenta and constrained by fits to the 
Near Detector f^-CC spectra measured at different 
beam configurations In contrast, the t'^^-CC spec- 
trum is less sensitive to the beam configuration since 
the majority of Dfj, arises from low-px parents (mainly 
7r~). Therefore, the pt shape of 7r~ parents is con- 
strained by the NA49 measurement of the relative 
7r+/7r~ production while the parent Pz shape and 
the absolute normalization are obtained from the fit to 
the Near Detector P^-CC spectrum. The result of the 
tuning on the antineutrino spectrum is demonstrated 
in Fig. m 



Figure 3: The three variables used for selection: the 
NC and CC selector (top), significance of the fit curva- 
ture (middle) , and the relative angle described in the text 
(bottom). Each variable is shown after selection on the 
others are applied. The vertical lines with arrow denote 
the optimized selection values. 



4. Results 

Fig. |5| shows the energy distribution of observed 
and expected I'^-CC events in the MINOS Far De- 
tector. We observe a total of 42 candidate events 
while expecting 64.6 ± 8.0(stat.) ± 3.9(syst.) events 
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Figure 4: The observed i>^-CC spectrum in the Near De- 
tector and the simulated spectrum before (blue) and after 
(red) tuning the meson production from the target (top); 
and the ratio of the predictions to the data (bottom). 



with no oscillation and 58.3 ± 7.6(stat.) ± 3.6(syst.) 
events with CPT-conserving oscillations when the an- 
tineutrino oscillation parameters are assumed to be 
equal to the best-fit MINOS neutrino oscillation pa- 
rameters (sin^(2^) = sin^(20) = 1 and Am^ — Am^ — 
2.43 X 10~^ eV^). The shaded blue band around the 
CPT-conserving prediction represents the systematic 
uncertainty, which is dominated by the relative nor- 
malization uncertainty between the two detectors, the 
muon momentum measurement from curvature, and 
the uncertainty in the contribution of downstream 
(decay pipe) production to the 9^ flux. The total sys- 
tematic uncertainty is less than 10% over the whole 
energy region. 

Extensive consistency checks were performed to 
make sure the observed deficit between the data and 
the prediction is not due to detector, reconstruction 
or selection inefficiency. The whole analysis procedure 
was performed with an alternative event selection us- 
ing independent variables and two different near-to- 
far extrapolation methods which resulted in a similar 
deficit and result. The track finding efficiency, in par- 
ticular for exiting tracks, was checked using stopped 
and through going cosmic muons. The muon charge- 
sign assignment was checked using an alternative track 
fitter. In addition, all far detector events with a /i+ 
candidate, events with a /i~ candidate ending close to 
the detector edge, and events with no reconstructed 
track were visually scanned for any sign of reconstruc- 
tion pathology. 

The effect of neutrino oscillation or transition can 
be applied to the Far Detector prediction using the 
parametrization of Eq. [T] or [21 respectively, and then 
fitted to the data to maximize the log-likelihood ratio 
with respect to the oscillation parameters or a. Confi- 
dence limits on the fit parameters are extracted using 
the Feldman-Cousins method (lo| with systematic un- 
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Figure 5: The observed j/^-CC spectrum (dots with error 
bard) and predictions in the MINOS Far Detector. The 
red histogram represents the prediction with no oscillation 
while the dashed blue with CPT-conserving oscillation. 
The shaded histogram is the expected background in the 
latter prediction. 



certainties incorporated. 

Fig. [5] shows the resulting contours on the muon an- 
tineutrino oscillation parameters, sin^(2^) and Am^, 
together with limits from a recent global fit to earlier 
data [ll| as well as MINOS limits on the correspond- 
ing neutrino oscillation parameters [l|. Although the 
best fit point is at high Ato^ because of the 1.9(7 over- 
all deficit, the data are consistent with CPT-invariant 
oscillation parameters at 90% C.L. and exclude the 
no-oscillation scenario at 99% C.L. We also perform 
a 1-parameter fit at maximal mixing (sin2(20) = 1), 
which helps to rule out previously allowed regions of 
Ato^ (Fig. [7]). At maximal mixing, /S.fh? < 2.0 x 10~^ 
and 5.1 x 10"^ < Am^ < 81 x 10^^ are excluded at 
90% conficence. 

There is no evidence for excess of events at low en- 
ergy that would indicate i'^ — > transition, and a 
fit to the data under this assumption gives an upper 
limit a < 2.6% at 90% C.L. 
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5. Future prospects 

MINOS has already collected about 7.2 x 10 
data in the i/^-dominated beam configuration, 
analysis of this data is under way and updated results 
are expected in the near future. 

More importantly for this analysis, the current of 
the NuMI focusing horns will be reversed in the Fall 
of 2009 in order to produce a dedicated antineutrino 
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Figure 6: Allowed oscillation-parameter regions from fit 
to the i/fi data (solid contours) together with limits from a 
global fit to prior data (dashed) and from the MINOS 
analysis (shaded). 



beam'^. Fig.[8]shows the expected D^j, and f^-CC spec- 
tra in the antineutrino mode relative the current neu- 
trino mode. The D^-CC event rate per PoT in the 
antineutrino mode is about a factor of three lower 
than the z^^-CC event rate in the neutrino mode due 
to lower production rate of negative mesons in the 
target and about a factor of two lower i?^-CC cross 
section in the detector. Despite this, the new configu- 
ration will significantly increase the i?^-CC event rate, 
in particular in the low energy region which is more 
sensitive to the oscillation, compared to the current 
beam configuration. 

Although the run plan is not yet finalized, we ex- 
pect to collect about 2 x lO^*' PoT data in the antineu- 
trino configuration in one year of running before the 
Summer shutdown in 2010^. Fig. [5] shows how much 
improvement can be achieved in sensitivity to the an- 



^The polarity of the magnetic horns was reversed and an- 
tineutrino running commenced on September 29, 2009. 

"^At the time of writing, the plan is to run in antineutrino 
configuration until March 1, 2010. 
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Figure 7: Confidence level of one-parameter fit to the MI- 
NOS data (blue) compared to that from the global fit 
(black). 
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tineutrino oscillation parameters with 2 x lO'^" PoT 
data collected in antineutrino mode compared to the 
existing 7.2 x 10^° PoT data in neutrino mode. This 
will allow us to measure the oscillation parameters, 
in particular Ato^, for muon antineutrinos with high 
precision and establish antineutrino oscillation at 5cr 
confidence level. 
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